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We describe a family of phase transitions connecting phases of differing non-trivial topological 
order by explicitly constructing Hamiltonians of the Levin- Wen[PRB 71, 045110] type which can be 
tuned between two solvable points, each of which realizes a different topologically ordered phase. 
We show that the low-energy degrees of freedom near the phase transition can be mapped onto 
those of a Potts model, and we discuss the stability of the resulting phase diagram to small pertur- 
bations about the model. We further explain how the excitations in the condensed phase are formed 
from those in the original topological theory, some of which are split into multiple components by 
condensation, and we discuss the implications of our results for understanding the nature of general 
achiral topological phases in 2 + 1 dimensions in terms of doubled Chern-Simons theories. 



I. INTRODUCTION 

The study of non-abelian phases of matter has 
drawn increasing attention in recent years, inspired in 
part by potential applications to topological quantum 
computatiorPEl. Since Moore and Read's^ proposal of 
the Pfafnan wave-function for the fractional quantum 
Hall state at v = 5/2, there have been many propo- 
sitions for realizing non-abelian matter in a wide va- 
riety of physical systems. These include in Sr^RuC^ 
superco nduct ors^^, Josephson junction arrayd^ES He- 
lium 3 A ! 14 * 15 !, cold-atom systems^^, as well as con- 
ventional superconductor interfaces on 3D topologi- 
cal insul ators^ and other strongly spin-orbit coupled 
material d 20 * 21 !. In addition to this wealth of new di- 
rections, recent experimental investigations of the 5/2 
state22H2U gj ve renewed incentive to study non-abelian 
matter. Despite this boom of interest, one topic which 
received relatively little attention until recentljEsHUl i s 
the question of phase transitions in these topological sys- 
tems, which will be the focus of this work. 

The non-abelian phases of interest can be character- 
ized by their topological order — that is, by the fusion 
and braiding (statistical) properties of their low-lying ex- 
citations. An interesting general question is to under- 
stand, on broader grounds, the possible transitions be- 
tween phases with different topological order. Drawing 
on the analogy with the Landau theory of symmetry- 
breaking phase transitions, we consider what features of 
such transitions, and the relationship between the phases 
they connect, can be deduced from the topological order 
alone. There are two questions to address here. First, if 
the phase transition is second-order, one would like to be 
able to deduce from the topological orders of the initial 
and final phases a long- wavelength description of the crit- 
ical theory. While we will not provide a complete answer 
to this question here, we will identify the critical theories 
for a large class of phase transitions, and comment on 
their universality. Second, we would like to understand 



- using information about the topological order alone - 
how the excitations of the phases on both sides of the 
transition are related. A substantial step in this direction 
was made by Bais and Slingerlancp2H3iJ who developed a 
framework known as topological symmetry breaking which 
describes topological phase transitions as a type of Bose 
condensation and deduces the topological order of the 
condensed phase. Here we will elaborate on this picture, 
by studying the fate of the excitations of a Hamiltonian 
which can be tuned through the phase transition. 

Our method for addressing the above issues will be 
to study a class of topological symmetry-breaking phase 
transitions which we can realize explicitly by a simple 
deformation of lattice Hamiltonians of the Levin- WerPS 
type. The key advantage of this approach is that, be- 
cause these Hamiltonians can be solved exactly at cer- 
tain points in the phase diagram, we will be able to 
understand in detail the phase transition (which in this 
case is always of the transverse-field Potts type) and the 
fate of the excitations on both sides of the phase bound- 
ary. Further, many of the features of these transitions 
are associated with the topological orders of the phases, 
rather than to the specific lattice model. Indeed, from 
the point of view of studying the general characteristics 
of the phase diagram, the Levin- Wen models should be 
viewed not as candidates for accurately describing a the- 
ory on the microscopic scale, but as "a-models" for these 
types of topological phases: though they will differ mi- 
croscopically from any known candidate physical system, 
they represent faithfully - in a highly tractable form- 
trie long-wavelength dynamics of any gapped phase with 
the same topological order. Hence by studying phase 
transitions in these models, we may hope to extract the 
generic features of the low-energy theory: the topological 
orders of the two phases, and (in some cases) the critical 
theory (or possible critical theories) of the phase transi- 
tion separating them. This is exactly as in the Landau 
symmetry-breaking program. 

This work is divided into three main parts. In Sect. [TTJ 
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we will show explicitly how to construct a lattice Hamil- 
tonian which can be tuned between two phases with dif- 
ferent topological order. The transition occurs via con- 
densation of a bosonic plaquette excitation which can 
be mapped onto a Potts spin; we will show that, un- 
der certain assumptions, the critical theory is that of the 
Transverse-field Potts model (TFPM). Sect. Ill discusses 
the fate of these results when the model is perturbed. In 
Sect. IV we will study the condensed phase in detail, 



showing how its excitations arise from those of the un- 
condcnsed model. O f prim ary interest here is the phe- 
nomenon of "splittin g ^ 30 ! 31 *' (in the context of conformal 
field theory this is related to fixed point resolution^) a 
single species of excitation into two distinct particle types 
after condensation, which we can understand explicitly in 
the lattice model. We conclude with a summary of our 
results and a brief discussion of other types of transitions 
which can be studied by the same methodology 



II. CONDENSING ACHIRAL SIMPLE 
CURRENTS 

In the present work we will realize phase transitions 
in the TFPM universality class by condensing a partic- 
ular type of boson (an achiral simple current) which has 
the property that multiple bosons combine (or fuse) ac- 
cording to addition rules appropriate to a Q-state Potts 
spin (a Zq spin). To describe precisely the nature of 
these transitions, we will study a family of lattice Hamil- 
tonians that are equivalent to exactly solvable Levin- 
WerPSl models at two points in the phase diagram. In 
the uncondensed phase, the topological order will be that 
of a "doubled" Chern-Simons theory^- a chiral Chern- 
Simons theory together with its mirror image. Here we 
will show that when the condensed boson $ is a certain 
type of achiral simple current with Zq symmetry (i.e., 
= Id), the long- wavelength description of the criti- 
cal theory can be mapped exactly onto the Q-state Potts 
model. (We will discuss in section III various pertur- 
bations to this Potts model). Condensing achiral cur- 
rents necessarily leads to a net achiral condensed phase, 
whose topological order is also captured by a Levin- 
WerPS model; the study of this phase will be the focus of 
Sect. ITvl 



A. Topological Lattice Hamiltonians 

Our starting point will be an exactly solvable Hamil- 
tonian of the type introduced by Levin and WerP^. We 
will restrict our discussion to models which realize dou- 
bled Chern-Simons theories, though the construction of 
Ref. |36] is more general. 

The Levin- Wen models can be viewed as deformed ver- 
sions of a lattice Yang-Mills theory. Their Hilbert space 
consists of a finite set of possible states i € {Id, . . . , r} on 
each edge of a honeycomb lattice. The labels {Id, . . . , r} 



of these states are analogous to the set of possible elec- 
tric fluxes in the gauge theory. (The label Id will always 
denote the trivial flux). The Levin- Wen Hamiltonian is 
constructed from two sets of commuting projectors: 
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The first sum contains projectors Vy acting on vertices 
V, and ensures that the ground states obey the con- 
straints V • E = at each vertex. For example, if the 
gauge group is SU(2), in the lattice Yang-Mills theory 
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where the rules for addition of angular momenta stipu- 



late that i X j = Y?i= J \i-j\ The projectors Vp^ in the 
second sum impose a condition analogous to V x B = 
at each plaquette P on the ground states. (As in lat- 
tice Yang-Mills theories, the operator Pp is essentially 
a superposition of Wilson loop operators encircling the 
plaquette P; we will will return to its the precise form 
presently.) All of the projectors in Eq. commute, since 
these two constraints can be satisfied simultaneously. 

The models are 'deformed' in the sense that the num- 
ber of fields r is finite, even though the gauge group is 
not discrete. For example, if the gauge group is SU(2) 
the lattice Yang-Mills theory would have electric fluxes 
corresponding to all allowed values 0, 1/2, ... of the total 
spin - of which there are infinitely many. In the Chern- 
Simons theory there is a maximum spin fc/2; the rules 
for adding angular momenta which specify i"V(or more 
generally, combining fluxes carrying different representa- 
tions of the gauge group) must also be modified to be 
consistent with this truncation, (k is referred to as the 
level, and the resulting Chern-Simons theory is denoted 
SU(2)k). Tabulations of these fusion rules for a number 
of theories can be found in Refs. 15514101 

It is important here that there is no U E 2 " term in the 
Levin- Wen Hamiltonian: there is an energy cost to cre- 
ating matter sources, but once a pair of these is created, 
there is no Coulomb-like energy cost associated with sep- 
arating them. This ensures both that the Hamiltonian 
is solvable, as all terms commute, and that the ground 
states and low-lying excitations are independent of all 
length scales, as required if they are to be topological. 

Because all of the projectors in Eq. ([I]) commute, the 
topological order of the model can be deduced from the 
statistical properties of its excitations. Though the elec- 
tric fluxes used to label the edges in the lattice are drawn 
from the representations of a chiral Chern-Simons the- 
ory, the Levin- Wen models realize a phase whose topo- 
logical order is that of a doubled (net achiral) Chern- 
Simons theory, whose sources we will label ip x with 
i,j G {Id, ..r} (indicating a composite excitation com- 
prised of a right handed i particle and a left handed j 
particle). 
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Excitations in the Levin- Wen model can be described 
in terms of electric sources, which violate the vertex con- 
dition Py — 1, and magnetic sources which violate the 
plaquette condition P p ^ = 1. As described in Ref. 2TJ 
these are related to the sources in the doubled Chern- 
Simons theory as follows. The electric excitations have 
the chirality of the original model (assumed to be right 
handed), so that electric sources are of the form i^x Idh 
or just in for short. The magnetic sources are achiral 
composites x of a particle and its mirror image. The 
"left-handed" excitations ij, carry both electric charge 
and magnetic flux of the particle type i. We will dis- 
cuss these excitations in more detail when we consider 
the topological order of the condensed phase in Sect. |IV| 

Importantly, while generically electric sources have 
nontrivial braiding properties, magnetic sources are al- 
ways bosons in the sense that they have trivial braid- 
ing statistics with themselves, hence can be condensed^. 
(In doubled Chern-Simons theories achiral particles are 
necessarily bosonic - meaning that they have trivial self- 
braiding statistics (in at least one fusion channel^}) 

Here we will condense a particularly simple type of 
magnetic excitation - one which behaves essentially like 
the magnetic flux in a discrete Abelian gauge theory. 
Specifically, we choose a particle <fi from the chiral theory 
with the property that <jfi = Id, <fi p ^ Id for any p < Q. 
Excitations with this property are known as Zq simple 
currents. In this work we consider condensation of the 
magnetic excitation $ = (f>^ x <fr L which is a bosonic Zq 
simple current carrying flux 4>. (Throughout this text, we 
will use <j> t° denote the chiral excitation, while <& refers 
to the corresponding achiral magnetic particle.) 

The particle we will condense is an excitation which 
violates only the plaquette projectors of Eq. (!]). To see 
how to condense such a particle, we must understand in 
more detail the form of the plaquette projectors. These 
have a very similar form to the plaquette term in a lat- 
tice gauge theory: they can basically be viewed as Wil- 
son loop operators around a single plaquette, which in 
lattice gauge theory contribute a "B 2 " term to the long- 
wavelength Hamiltonian. Since we wish to work in a basis 
where edge labels denote electric flux, the result is an op- 
erator which raises or lowers the electric flux on each edge 
of the plaquette. This will be familiar to many readers 
from the standard formulation of Ising gauge theory, in 
which the magnetic term in the Hamiltonian is commonly 
expressed as an operator which simultaneously flips all 
spins bordering a plaquette. 

The amount by which the electric flux changes under 
these operators is determined by the representation car- 
ried by the Wilson line. To construct a plaquette projec- 
tor P") onto magnetic flux j, we must use a particular 
linear combination of these representations, given by: 



(3) 



i=id 



in the language of Ref. |36J) around the plaquette P in 
the representation i, and the sum runs over the (finitely 
many) representations i. Here S is the modular S matrix 
of the chiral Chern-Simons theory and is related to the 
braiding properties of the excitations therein. The Levin- 
Wen Hamiltonian (Eq. [l| has projectors Vp^ on each 
plaquette which project onto the identity representation 
Id, in which case 
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(4) 



where Aj is a real number known as the quantum dimen- 
sion of the representation i, and T> = Jd A?. 

More generally, if is a Zq simple current, then 5^ t 
has the particularly simple form 



(5) 



where % is specific to <f>, but independent of n. (To be 
precise, qi is an integer, with qi = if i and <j> have triv- 
ial mutual braiding statistics in the chiral Chern-Simons 
theory). Substituting |5]) into ([3| gives a projector onto 
plaqucttcs with flux (j). 



B. Condensation 

We now describe the simple deformation of the ex- 
actly solvable Hamiltonians which allows us to tune the 
system through a phase transition in which the boson 
$ = <f>n X 4>l condenses, and derive a second exactly 
solvable Hamiltonian which captures the topological or- 
der of the condensed phase. Here we restrict ourselves to 
the case that $ is a purely magnetic excitation, and that 
it is a simple current. 

We first require an operator that will pair-create the 
requisite vortices. If is a Zq simple current, we may 
create a pair of vortices of flux <fi n and <t$~ n on adjacent 
plaquettes P± and Pi by acting with an operator: 



vUr)\h2) 
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(6) 



where Wi(P) denotes a Wilson line (or "string operator" 



where ei2 is the edge between plaquettes Pi and P2, and 
|i 12 ) denotes any state that has electric flux i on the 
edge ei2 shared by plaquettes Pi and Pa. (Which pla- 
quette obtains a vortex of flux </>", and which a flux <ffi~ n , 
depends on which orientation we choose for the electric 
flux on the edge i.) This operator 2 (4> n ) thus assigns 
a phase to the wave-function depending on the value of 
the electric flux on the edge separating plaquettes Pi 
and P2 (Fig. [T]) . Note that if the plaquettes already have 
magnetic fluxes <fi a and cj) b , then this operator simply in- 
crements these fluxes mod Q accordingly to ^(™+ a ) mod Q 
and ^(<2-™+ b ) mod Q. 

To engender condensation, we will simultaneously 
decrease the energy gap of excitations (magnetic 
fluxes of 4> n ) without changing the gap to the model's 
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FIG. 1: Creating a pair of simple-current vortices on ad- 
jacent plaquettes. The operator Vj 12 ((p n ) acts on the edge 

separating plaquettes Pi and P%, assigning a phase e 2 ^ 1 ^ 
to any component of the wave-function in which this edge is 
labeled by i. If the plaquettes start without flux as shown in 
the figure, the result is a pair of vortices of flux <j> n on Pi and 
(j> Q - n on P 2 . 



other excitations, and also add a term proportional to 
V'(<p) which creates such excitations spontaneously. The 
Hamiltonian that does this has the form: 



(7) 
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E 



V 



(0) 



P n=0 



p 



(Note that the sums over powers of (f> in the first and 
second lines run from to Q — 1 and 1 to Q — 1, respec- 
tively). Here is the edge shared by plaquettes Pi and 
Pj, and 



1 



(8) 



n=0 



_1 with equal amplitude. 



creates all fluxes 0, <f>, <f> , ... , 
Note that here V} = V e .. = 

e z3 ^ji 

In the limit ax — l,ctjv — 0, this Hamiltonian is 
just the exactly solvable Levin- Wen Hamiltonian (Eq.[T]), 
which assigns a mass gap of e e to electric sources, and e m 
to magnetic sources. The second line of Eq. [7] allows us 
to tune the model through a condensation, by varying 
the parameters oln and ax- 

To understand the phase portrait of ([7]), it is helpful to 
consider the effect of the various terms on the subspace 
of states in which the only excitations are the magnetic 
<I> excitations which we wish to condense. The first line 
of (JtJ) is indifferent to the presence of 4>, as $ does not 
violate the vertex condition, and the sum of plaquette 
projectors 



Q-i 
E^ 



(<T) 



(9) 



has eigenvalue 1 for a plaquette excitation $ n (magnetic 
flux (f> n , for n = 0, 1, ...Q — 1), and for a plaquette con- 
taining any other possible excitation. The role of these 
terms of the Hamiltonian is to ensure that there is a min- 
imum gap of %p for all magnetic excitations other than 



and ev for all electric excitations, everywhere in the 
phase diagram. (Here we keep e m and ey fixed). These 
energy gaps are unfrustrated, in the sense that neither 
the plaquette term nor the creation term on the second 
line have any amplitude to create these other excitations, 
which therefore remain both gapped and orthogonal to 
the ground state throughout the phase transition we de- 
scribe. Thus we will be able to focus our entire attention 
on the Hilbert space containing only by $ n type excita- 
tions. 

The second line of Q contains the terms which drive 
the system through a condensation transition. The first 
term 



~2~ aT E 



,(o) 



Q-i 



V 



(10) 



sets the scale of the gap to creating a $ excitations, which 
we may tune from e m at ax = 1 to at ax — 0. For 
clt < the formation of $ excitations is energetically 
favored; for small negative ax the system remains in 
the condensed topological phase. If, however, we make 
ax < 0, \ax\ 3> oln we will find a new variant of the un- 
condensed phase in which each plaquette is occupied by 
a definite superposition of <fi n vortices. (Taking a at = 
in this phase gives another exactly solvable model, again 
of the general form discussed in Ref. 36). 

The term in ([7| with coefficient ajv adds an amplitude 
to create or destroy pairs of fluxes on adjacent plaquettes 
(as shown in Fig.[T]). Analogous to adding an anomalous 
term to a Hamiltonian of conventional bosons 



H = H a + A* b ] + Ab, 



(11) 



spontaneous creation of $ particles is precisely what we 
should expect to need to add to a Hamiltonian in order to 
form a $ condensate. In these topological models, par- 
ticles must be created in particle-antiparticle pairs con- 
nected by a Wilson line (or string operator) ; Vj. . acts on 
pairs of neighboring plaquettes to create these. 

It is instructive to consider the Hamiltonian ^ in the 
limit ax — 0,ajv = 1. The plaquette term then has the 
form: 



O-i 

n=0 



r /O-l \ 



= E <WA^(P) 



(12) 



j=ld 



where Wj(P) is the Wilson line (or string operator) 
around plaquette P carrying the representation j. For 
representations which have trivial braiding statistics with 
the 4> flux, qj — 0, and the second sum on the right hand 
side (in the parentheses) just gives an overall factor of 
Q. If qj ^ (mod Q), then the second sum leads to 
complete destructive interference, and the Wilson line j 
is eliminated from the plaquette term. 
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Similarly, the effect of the vortex creation term V on 
the edge label i on edge e is: 

, O-i , O-i 

nE^l^nE^I^W') • ( 13 ) 



Hence — ajve m V e effectively assigns an energy cost ajve m 
to any (electric flux) label i on edge e for which 7^ 0. 
This creates a confining potential for any label which 
braids non-trivially with the flux (j) of the condensed bo- 
son $ = (j) R X 

In addition to the Levin- Wen point ay = 1, oln = 0, 
we therefore have a second special point in the phase dia- 
gram at ax = 0, oin > where the Hamiltonian is again 
exactly solvable. This point represents the exactly solv- 
able (and fully topological) limit of the condensed phase, 
in which labels with qi ^ (which do not braid trivially 
with <f>) have been completely eliminated from the theory. 
This can be consistently done because the plaquette term 
(which now contains only raising operators with q t = 0) 
does not mix edge labels with different qi. Hence states 
containing only edge labels with q t = comprise the en- 
tire low energy space of the Hilbert space in this limit. 

Restricted to these states, at olt = 0,a!Ar = 1, the 
Hamiltonian |7]) has a particularly simple form. Since 
V e = 1 on all remaining states, we may drop it from the 
Hamiltonian, leaving 



H eff = 



(0) 

p 



where 



AiWi(P) 



(14) 



(15) 



has the general form of a plaquette projector onto flux, 
as per Eq. M. As Vp^ and Vy commute, Eq. (14 1 is 
again a Hamiltonian comprised of commuting projectors, 
of the same general form as our initial Levin- Wen Hamil- 
tonian. The difference here is that the labels i for the 
edges are now drawn from the subset of the original la- 
bels for which = 0. We will study the consequences of 
this restriction on the model's excitations in Sect IIVI 



C. Effective Potts model of the phase transition 

We have thus argued that (JtJ) can be tuned through 
a phase transition between two solvable Hamiltonians of 
the Levin- Wen form, with each of the solvable models 
capturing exactly the topological features of one of the 
two phases. (See Ref. 02] for a discussion of an analogous 
family of achiral solvable lattice Hamiltonians inspired 
by Kitaev's toric code models^). Here we will exploit a 
natural mapping between the long-wavelength dynamics 
of Q and the Potts model to identify the phase transition 
as that of the transverse-field Potts model. 



To carry out this program, we define an effective low- 
energy Hilbert space which we will map exactly onto the 
Hilbert space of the Potts model. Specifically, since all 
excitations other than <J> n remain gapped (with a mini- 
mum energy of min{ey, %*-}), we need only consider the 
subspace of states containing just the ground states of the 
initial Levin- Wen model, and states which can be derived 
from these by adding some number of magnetic $ (and 
$ n ) excitations. We emphasize that this subset is closed 
under the dynamics of Eq. Q for all parameter values 
of the Hamiltonian, and this subspace contains all modes 
which become critical at the phase transition. Hence this 
low-energy subspace can be consistently separated from 
the rest of the Hilbert space in order to understand the 
dynamics of the critical point. 

Because there is a unique way to combine such exci- 
tations ($ a x $ b = $(o+f-)m°d<3) 5 a state in the Hilbert 
space is uniquely determined by specifying the topolog- 
ical ground state sector and the number np of <fi fluxes 
on each plaquette P. Hence we may define the basis 



\g;ni,n 2 , ...n N ) 



(16) 



where and € 0, . . . , Q — 1 denotes the number of <f> flux 
quanta through plaquette i (i.e, the plaquette contains 
excitation <f> ni ). Here, g labels a ground state sector of 
the unperturbed Levin- Wen model (Eq. [I]) in cases where 
the manifold has handles. Our mapping to this Q-state 
model on plaquettes (i.e., on the dual lattice) is illus- 
trated in Fig. [2] a. 

Readers familiar with the general theory of TQFT's 
should note that a state of the form ( 16 ) uniquely iden- 



tifies a state in the original lattice model because $ is 
a simple current. For other types of magnetic excita- 
tions, such as Fibonacci anyona^l, to specify the state of 
the lattice model uniquely requires additional informa- 
tion about the relative fusion channels of the excitations. 

The terms in the first line of ^ act as the identity 
on our effective Hilbert space. In the basis (16), the two 



terms in the second line have the matrix elements: 

(nWvfl - Pjf Vi) = <k,n a (2<W - 1) (17) 

71=1 

(ni,r4|V r e t 13 (^ fc )|ni,n 2 ) = S n > i± ^ ni 5 n ^ n2 (18) 

with the ± in the second line given by the chosen ori- 
entation of the bond e\2- In this reduced Hilbert space, 
below the gap to excitations out of this Hilbert space our 
Hamiltonian Eq. [7] is then precisely 

c 



spin 



- a T^^(2<5„ l , -l)|n. i )(r 



(19) 



E E i n « + k - n j ~ k )( n i> n ]\ 



<i,j>rii,nj k 



The Hamiltonian Q keeps the gap to all powers n = 
1, ... Q — 1 of $ the same, reflected in the first term in 
_ff S p m . The second term gives equal amplitudes for tran- 
sitions that increment and decrement fluxes by any k, in 



accordance with the creation operator Q . These choices 
will give us precisely an effective Q-st at e transverse-field 
Potts model (TFPM) of the phase transition. In section 
|III A| below we will consider how we may change mod- 
ify -Hgp m and obtain many possible Hamiltonians for Zq 
Potts-like spins (for example, the Q-state clock model), 
and discuss how this may affect the phase diagram. 

To obtain a more convenient representation of the 
TFPM, we re-express (HTfb, (fl8|) in the basis 



Q-i 



(20) 



n=0 



The matrix elements ( 18 ) become precisely the Potts in- 
teraction: 

(liX\Wjli,h) = 8 Vlth 6u l , h 5 ll - la ,o (21) 

while ( fi"7| ) becomes the appropriate transverse-field 
ternPl 



O-i 



(22) 



This maps the long-wavelength description of the phase 
transition precisely onto the conventional representation 
of the transverse-field Potts Hamiltonian: 

tfspm = ~ «tE fym-mui) 

-e m aN fik,lj\k,lj}(k,lj\ • (23) 
<ij> 

(Our choice of notation now becomes clear: is the 
transverse field term and a at is the neighbor interaction) . 
Here ajy > is used to promote condensation, so the 
model has a ferromagnetic phase in which all of the l- 
spins are aligned and a paramagnetic phase where they 
align with the transverse field. 



In terms of the Potts model description, the labels 
which become confined in the condensed phase corre- 
spond to domain walls of the ferromagnetic Potts model. 
Specifically, if the edge separating plaquettes i and j 
is labeled in, then the Potts interaction is 



(24) 



Hence if g ei ^ 0, then in the basis of Eq. ( 20 ) , Zj — lj ^ 



and there is a domain wall in the ferromagnetic Potts 
representation. Thus these labels become confined since 
domain walls are confined in the ferromagnetic phase (See 
Fig. [2]b). (Though the full Hilbert space of the lattice 
model also includes open strings, for ey > these do not 
occur in the zero-temperature phase transition, due to 
the fact that all terms in H commute with P\r)- 




o co oo 



FIG. 2: Effective mapping from the low-energy sector of the 
string-net model to the Potts model, illustrated here for the 
Ising case, (a) In the low-energy sector, we retain only mag- 
netic excitations of flux <f> n , as all other excitations remain 
gapped throughout the phase transition. Up to an index spec- 
ifying the topological ground state sector, the relevant states 
are specified uniquely by assigning a Potts spin n £ {0...Q — 1} 
to each plaquette. In this figure, an arrow pointing right 
(left) indicates n = (1) on that plaquette. (In the spin 
model we identify these with S x = 1 and —1 respectively.) 
(b) The terms in the second line of the Hamiltonian Eq. Q 
act non-trivially on these states: the plaquette term controls 
the energetic cost of creating a flux, which we identify with 
the Potts transverse field. The flux-creation term V e gives a 
ferromagnetic Potts interaction in the basis of Eq. ( |20[ ), indi- 
cated here by arrows pointing up (down) to denote the states 
I = (I = 1) on that plaquette. (or S% = 1 and —1). The 
eigenvalue of V e is 1 if the edge e is labeled by a represen- 
tation i with qi — 0, and otherwise - indicating that any 
edge label i with qi ^ (mod Q) signals a domain wall in the 
Potts model. 



There is one important difference between the lattice 
model and its spin analogue: while the latter has Q 
symmetry-related ferromagnetic ground states, the for- 
mer has only one. The topological model is only sensitive 
to the locations of the domain walls in the Potts model, 



not to the orientation of the Potts spins in the basis ( 20 1 



(This loss of ground-state degeneracy also occurs in the 
dual loop- gas representation of the Ising model). This is 
because the Hilbert space of ^ has the restriction that 



rii = mod Q. 



(25) 



That is, the Hilbert space we use is not identical to that 
of the Potts model; rather, since vortices can only be cre- 
ated in vortex- ant i- vortex pairs in the topological theory, 
it contains only states connected to the vacuum by the 
action of the Potts interaction (or pair-wise spin flips, 
in the Ising case). Since the spin Hamiltonian connects 
only states within this subspace, this difference does not 
affect the system's dynamics or the nature of the phase 
transition. However, the Q-fold the degeneracy of ferro- 
magnetic ground states in the Potts model is absent in 
the topological model: in the latter, there is no operator 
which measures the eigenvalue of I on a particular pla- 
quette, since such an operator would necessarily change 
the number of vortices on only one plaquette. We note 
that this difference renders certain interesting modifica- 
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tions of these models^^ difficult to realize in the present 
context. 

We may therefore conclude that the phase transition 
of ([7]) is precisely that of a ferromagnetic transverse-field 
Potts model in 2 + f dimensions, since we have con- 
structed an explicit mapping between the two. The orig- 
inal string-net model maps to the paramagnetic phase 
olt S> Qijv, where the transverse field dominates to 
ensure that vortices are relatively rare. In the ferro- 
magnetic phase of the Potts model, the vortex creation 
term dominates in Eq. ([T]). In this regime <&" vortices 
have condensed and certain edge labels become confined. 
The transition between the two will be in the 2+1D 
transverse- fiel d-Po tts universality class (which is first or- 
der for Q > jPEZI). 

The notion that a phase transition which changes the 
topological order can be described by a spin model, in 
which a global symmetry is broken when a local vari- 
able obtains an expectation value, is somewhat counter- 
intuitive. To understand this correspondence, notice first 
that after projecting onto the states relevant to the criti- 
cal theory, we arrive at an effective description in a pure 
(matter-free) discrete Abelian gauge theory. That is, the 
vortices that we condense behave exactly like vortices of a 
discrete Abelian gauge theory, both in terms of how mul- 
tiple vortices combine, and from the way in which the 
operators which create them act on the gauge-invariant 
states. The mapping to the spin model then simply ex- 
ploits the fact that this discrete Abelian gaug e theo ry 
without matter sources is dual to a spin model^^l ( a 
model of matter without magnetic vortices). This du- 
ality has been exploited previ ously to study phase tran- 
sitions in the Toric cod e! 50 * 51 -! Once matter sources are 
added to the gauge theory, the dual theory is again a 
discrete gauge theory: since itinerant charges will always 
feel the Berry's phase of the condensing vortices, we must 
include gauge fields in the dual spin model such that the 
true gauge- invariant order parameter is non-local in both 
representations. 



1 is ever obtained: 



D. Example 



To make the discussion of the previous section con- 
crete, let us give here an example. Other, more general 
examples are discussed in Sect. IV D 



First, as described in Ref. [HJ we may construct the 
initial Hamiltonian from representations of the Chern- 
Simons theory SU (2)2- In this case there are three pos- 
sible values 0(= Id), 1/2, and 1 of the total spin. The 
(commutative and associative) fusion rules of these spins 
are similar to angular momentum addition, except that 
the rules are truncated such that no value greater than 



Oxj 
1/2 x 1/2 
1/2 x 1 
1 x 1 



3 

+ 1 

1/2 





(26) 



Note that the particle 1 here is a Z2 simple current. 

The vertex projectors of the Levin- Wen model, based 
on these fusion rules, have eigenvalue 1 if the edges in- 
cident on a vertex are in one of the following 3 com- 
binations (here the order in which they appear is not 
important): 

(0,0,0) (0,1/2,1/2) (0,1,1) (1/2,1/2,1) , (27) 

and the vertex projector gives zero for any other combi- 
nation of edge labels. 

The Wilson lines (whose action is also based on the 
rules (26 )) raise and lower the edge labels according to: 



W \i) = |t) Wi\0) = |t) 
W 1/2 \l/2) cx a\0) + W 1/2 \l) cx |l/2> 

Wi|l/2}oc|l/2) Wi|l)oc|0) (28) 

Here the constants of proportionality depend in general 
on the labels of adjacent edges, as well as the edge being 
acted upon (See Ref. [551 for details); however, their pre- 
cise value is unimportant for our purposes. The last line 
reflects the fact that the spin 1 particle is an order 2 (Z2) 
simple current: raising the edge label 1 by W\ necessarily 
gives the trivial label 0. 

We condense the spin 1 magnetic excitation (IrX 1l), 
which is a Z2 bosonic simple current. As described in the 
previous section, the effective theory below the gap can 
be mapped precisely onto that of a 2d transverse-field 
Ising model. To do so, we assign a spin variable S x = — 1 
to every plaquette containing a spin-1 vortex, and S x — 1 
to vortex- free plaquettes. The relevant components of the 
S matrix have the form: 

5*1,0 = Si,i = 1 S Q ,i = \/2 S h i = -V2 . (29) 

In particular, we have: 



(30) 
(31) 



and the term which creates a pair of vortices (and hence 
flips the spins) on adjacent plaquettes is 



V eij = (-l) 2 ^ = SfS 



(32) 



where s measures the total spin of the representation on 
the edge e^ . Thus we identify the spin-1/2 label with 
the domain wall in the Ising model: an edge carrying the 



spin-1/2 representation necessarily separates two plaque- 
ttes with oppositely oriented Ising spins. 

In the solvable limit of the condensed phase, we are 
left with only the edge labels and 1. The vertex condi- 
tion is now that an even number of edges labeled 1 must 
enter each vertex, while the plaquette projector (301 flips 



the label on all edges surrounding the plaquette from 
to 1 or vice versa. Some readers may recognize this as 
the Hamiltonian for Kitaev's Toric code (or Ising gauge 
theory with matter). We will return to this point in Sect. 

El 



III. UNIVERSALITY IN TOPOLOGICAL 
SYMMETRY BREAKING TRANSITIONS 

Thus far, we have established a precise mapping be- 
tween a family of Hamiltonians with the special prop- 
erty that they are exactly solvable at two points in the 
phase diagram, and (at energies below the minimum gap 
to excitations at these solvable points) the transverse- 
field Potts model. However, our choice of Hamiltonian 
|7]) is quite non-generic: we would expect that any real 
physical system whose long-wavelength dynamics are de- 
scribed by the lattice model would include perturbations 
away from the solvable Levin- Wen limit everywhere in 
the phase diagram. Hence to apply our understanding 
of the phase diagram of the topological lattice model to 
more general systems, we must understand the effect of 
generic perturbations on its behavior. 

Ideally, we would like to achieve a framework analogous 
to the Landau approach to symmetry-breaking phase 
transitions. There are two important elements to this 
analogy. First, a phase can change its symmetry only by 
a phase transition. Specifically, in a gapped system small 
perturbations to the Hamiltonian which do not close the 
gap necessarily leave the symmetries intact. (The excep- 
tions to this rule are systems which undergo first order 
transitions, or second-order transitions associated with 
the formation of microscopic domains^ 1 -, in which case 
the gap itself need not close, but there are other sin- 
gularities associated with the phase boundary). Second, 
near a second-order phase transition the long- wavelength 
behavior is largely determined only by the symmetry be- 
ing broken. Specifically, symmetry dictates the nature of 
the field theory at the transition. There can be specific 
choices of parameters for which not all relevant opera- 
tors compatible with the symmetries of the two phases 
appear in the critical theory; however, these non-generic 
systems lie at unstable critical points which flow to dif- 
ferent critical theories if the perturbation is added. 

We emphasize that though the present work treats 
transitions in which a simple current (or Potts spin) 
condenses, these questions are relevant to more general 
TSB transitions. It is also possible to deform Levin- Wen 
Hamiltonians (by adding appropriate analogues of V e , 
and modifying the plaquette projectors appropriately) to 
condense other types of bosonic vortices, such as the Fib- 



bonaci anyons for which a transition of this type has been 
discussed by Ref. [32] This can be done without closing 
the gap to electric sources, so that again we may imag- 
ine studying an effective "spin model" which describes 
only the ground states and condensing vortices. (The 
resulting spin models will not be Potts-like, however). 
However, not all such constructions yield a second set 
of parameter values (analogous to or = 0, chn — 1) at 
which the Hamiltonian is again exactly solvable; hence 
it is more difficult to establish the properties (or exis- 
tence) of the condensed phase. To grasp completely the 
phase portrait of achiral topological phases would require 
an understanding of these more complex transitions and 
their universality, which we do not undertake here. 

The notion that a system cannot change its symme- 
try without undergoing a phase transition has a well- 
understood analogue for topologically ordered phases. 
Specifically, it is knowipSES] that any deformation to 
a topologically ordered phase which does not close the 
quasi-particle gap cannot change the topological order. 
Working within the Hilbert space of the model at hand, 
this means that we may include small perturbations 
which admix some number of electric or magnetic sources 
with the ground states, which nonetheless do not change 
the topological ground state degeneracy or long-range 
braiding statistics that we associate with the topologi- 
cal order of the solvable Hamiltonian. Hence topological 
order in the present context plays the role of symmetry 
in the more familiar Landau paradigm. (Since topolog- 
ical orders can often be associated with gauge theories, 
this is very natural. The symmetry broken at the phase 
transition, in these cases, is simply a gauge symmetry 
rather than a global one). 

The second question - of whether the critical theory 
is robust against such perturbations - is less well under- 
stood in the topological context. In the previous section 
we chose a special trajectory through the phase transi- 
tion in which the degrees of freedom implicated in the 
phase transition could be mapped exactly to those of 
a transverse-field Potts model. This suggested that the 
critical theory is dictated by the nature of the condensing 
boson. Specifically, when this obeys Zq fusion rules, we 
expect to find critical theories in the universality class of 
a Q-state spin model (a Potts model or similar). To un- 
derstand whether this conclusion about the critical the- 
ory is valid more generally, we must consider the effect 
of generic perturbations to the Hamiltonian ([7]). 

It is convenient to separate the possible perturbations 
into two classes. For Q > 3, it is possible to introduce a 
special type of perturbation which changes the parame- 
ters of the Potts Hamiltonian ( 23 1 . As we explain below, 



if Zq contains proper subgroups (i.e., for Q not prime), 
this results in a richer phase diagram than the one de- 
scribed above. Other types of perturbations may not 
have analogues in the Potts description; our chief con- 
cern here we will be to consider their effect on the critical 
theory perturbatively. 
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A. Perturbations within the Potts model subspace 

Let us begin by considering perturbations whose effect 
is to change the relative strengths of the ferromagnetic 
couplings or transverse-field terms for the various pow- 
ers of <j) in Eq. ((7|. For Q > 3 this allows considerable 
scope to modify the spin model, leading to a rich phase 
diagram. A thorough understanding of the behavior of 
these more general spin models is in itself an interesting 
question, whose general features we will outline here by 
reference to the existing literature where possible. 

We consider perturbations of the form 



SH = J2$H P -J2SH eij (33) 
P 

shp = eVh^ ^ = iE4 e V4(« 



n=0 



n=0 



where 

» _ 



c (P,e) 



can be chosen arbitrarily provided that: 1) 
er>-ra (required for hermiticity) ; 2) the sign of 



the net transverse field term for any power of <fi remains 
positive; and 3) that all Potts interactions remain ferro- 
magnetic. Provided these three conditions are met, the 
perturbed Hamiltonian will still undergo (one or more) 
TSB-type phase transitions. Here we will consider the 
case where e« and en are real, though the complex case 
has also been studied in the context of spin models^. 



Since the perturbations ( 33 ) clearly have no impact on 



the mapping to Q-state spins, we may consider the effect 
of these perturbations within the spin picture. In terms 
of the spin state \np) on the dual lattice (recall that a 
state | np) indicates a vortex particle $™ on the plaquette 
P in the original model), the operators in Eq. (33) are 
mapped to the operators: 



^spin = q E E 4 e) K-"> n J+ n )( n i> n il 

<i 7 j> 7ii ,rij ,n—0 



i rii—0 



(34) 



The first term flips pairs of spins on neighboring plaque- 
ttes; the second is the (transverse) magnetic field. To 
make the correspondence to the Potts and clock models, 



we re-express these in the basis (20), where the Hamilto- 
nian has the general form: 

(l[l' 2 \SH e \hl2) = 5 ViM 8 V2 ^e^co, 2m{ll Q h) 

n—1 

Q-i 

(l'\6H P \l) = Y e^e 2wi(, ' -,) t (35) 



The Q-state clock model is obtained for Cn' P ^ — 

S n i; the Potts model occurs for e„ = l,ei pS> = l,ei P ^ 



— l(n > 0), as described above. More generally, we may 

consider any nearest neighbor Zq spin interactions. 

Tuning the magnetic field leads to an even wider range 
(p) ( p) 

of possibilities. If e„ = £q_„ then the field is again a 

superposition of clock-model type fields in different di- 

(p) (p) 

rections. If e„ ^ e o- n the model is chiral, in the sense 
that the transverse field has either a left- or a right- ro- 
tating eigenstate (but not both). This can produce, for 
example, a special case of the chiral Potts model first de- 
scribed by Ref. [56] While we do not expect the choice 
of transverse field to affect the symmetry-broken phases, 
it could have an impact on the loci and nature of the 
phase transitions, which we will not attempt to describe. 
We note, however, that there is no such freedom in the 
Ising or Q = 3 Potts case - and that for Q > 2 the Potts 
transition is in any case first order in 3D. 

We therefore focus on the various possible interaction 
terms, and briefly describe the different possible phases 
allowed by symmetry-breaking in the spin model. There 
are two possibilities: if Q is prime, then condensing $™ 
for any n will result in ($ fe ) ^ for all k. In this case 
all ferromagnetic symmetry-broken ground states in the 
spin model are related by a global Zq rotation. The 
Potts Hamiltonian ^ is special since it represents the 
unique choice of transverse field term and ferromagnetic 
coupling corresponding to a Hamiltonian that is exactly 
solvable in the limit ax = 0, > 0. However, there 
are many choices of Hamiltonian which all undergo tran- 
sitions breaking the same Zq symmetry. Thus even for 
vanishing transverse field, there is a large family of pa- 
rameters for which the Hamiltonian is in the same 
phase but not exactly solvable. For Q > 3 the clock 
model is one such Hamiltonian. 

If Q is not prime, then there will be some choices of k 
for which we can have both ($ fc ) ^ and (<&) = 0, giving 
two or more distinct phases in which different subgroups 
of Zq have been broken. In the Potts basis of Eq. ( 35 1 , 
this means that 



(0 = «o 



Q/k 

E 

3=1 



(•■*) 



J I 



(36) 



where we have taken the simplest case = vq for all 

j = 1 . . . Q/k. Because Q/k is an integer, this results in 
a configuration of expectation values which is invariant 
under shifts in I by a subgroup of Zq. 

To illustrate this possibility, we take Q = 4. We may 
individually set the values of the transverse field and fer- 
romagnetic interaction for the two Z4 fields $, <I> 3 (whose 
couplings must be the same for the Hamiltonian to be 
hermitian), and the Z2 field <I> 2 . The ferromagnetic cou- 
pling of the spin Hamiltonian is 

h = j 2 j2 + + J 4 E ( 3? ) 

<ij> <ij> 

(Note that we are considering $ 3 , <J> 2 and $ to be three 
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independent fields here). Taking 



V2 



(S-icr) $ 3 



(38) 



we see that Eq. (37) is equivalent to the Ashkin- Teller 
model: 



H = J 2 y^ffjiSj + CiCTj) + J4 <7i<7jSiSj 



If we choose the transverse field term 



H' 



B 1 [s\ x) 



(39) 



(40) 



then the quantum problem in non-vanishing transverse 
field is equivalent to the 3D Ashkin- Teller model. 

The phase portrait of the classical 3D Ashkin- Teller 
model has been studied, for example, in Ref. |S7) If 
all couplings are ferromagnetic, there are two distinct 
phases: for J 2 7^ the full Z4 symmetry is broken, 
whereas If J 2 — 0, then there is a "<I> 2 " phase with 
(<I> 2 ) 7^ but ($) = 0. There are two Ising-like second or- 
der phase transitions separating the $ 2 phase from both 
the paramagnet and the fully symmetry-broken phase. 
The two phase boundaries end at a tricritical point, af- 
ter which the phase boundary separating the paramagnet 
from the fully broken phase is first order. 

In summary, there are a wide range of Hamiltonians 
which will ultimately condense the same simple current in 
the topological lattice model, corresponding to different 
choices of ei P ' e ' in Eq. (33 1. Equally, if Q is not prime we 
may condense simple currents of order Q' = Q/n, leading 
to new distinct symmetry-broken phases. 



B. 



Perturbations outside of the Potts model 
subspace 



Perturbations of the form ( 33 1 are special in that they 



do not introduce any excitations other than $ (Tl) vortices 
into the ground states of the system. In this case the task 
of understanding the phase transition reduces to one of 
understanding a spin model. This statement is, at sec- 
ond glance, a rather surprising one: we have reduced a 
question about topological orders, where the long-ranged 
statistical interactions mediated by gauge fields dictate 
the inter-particle interactions, to one about a spin model 
in which there are no gauge fields at all. Essentially this 
is because the magnetic particles we condense behave like 
vortices in an abelian gauge theory, so that we may ex- 
ploit an electric-magnetic duality to map the theory of 
vortices on the lattice onto a theory of 'charges' (here Zq 
spins) on the dual lattice. 

A generic perturbation to the Hamiltonian |7| will in- 
troduce excitations outside of this pure $™ vortex sector 
into the ground state, however. Here we consider the 



question of whether such perturbations qualitatively al- 
ter the phase diagram or critical behavior. We present 
general arguments that terms which weakly mix electric 
excitations with the ground states do not alter the critical 
behavior, although for Q > 2 the pure Potts transition 
is in any case first order. Perturbations generating other 
types of vortices correspond to annealed disorder in the 
(classical 3D analogue of the) spin system, which we also 
expect to be irrelevant at the critical point. 

We begin with a slightly perturbed version of our 
Hamiltonian: 



H = H 



LW 



H W + eH (E) 



(41) 



Here, Hlw is the original Levin- Wen model, and 
includes the terms which tune the system through the 
condensation transition (for example, those given in Eqs. 



IIIA). 



[7J or any perturbation discussed in Sect, 
contains perturbations which cannot be described in the 
effective spin model. The possible perturbations fall into 
three classes: H^ E ' may create deconfined electric exci- 
tations, magnetic excitations (which are always decon- 
fined), or confined electric excitations. We consider each 
possibility in turn, to determine whether any of these are 
relevant perturbations at the critical point. 

Sources of deconfined electric excitations have no effect 
on the Zq spin state; they act as the identity operator 
on the Zq spin model subspace. Hence a perturbation 
which excites only such sources merely changes the short- 
distance characteristics of the topological ground states 
over which the transition occurs. The altered ground 
states will still be compatible with our mapping to the 
spin model, however; hence the critical theory will be 
unaffected by this perturbation. 

Perturbations exciting magnetic sources v 7^ </>™ also 
fall into two categories, depending on whether (jf x v = v 
for some r < Q. These map onto dynamical dilutions in 
the transverse field and in the sites of the spin model, 
respectively. Specifically, suppose a plaquette P is occu- 
pied by a generic vortex with flux v in the uncondensed 
ground state. Acting with the operator which creates 
$ on this plaquette will produce a new vortex with flux 
v' — (f> x v. Since the product is unique, we may map this 
state onto a spin state by identifying v' on this plaquette 
with the Zq spin normally associated with <j>. The resul- 
tant mapping onto states of a spin model is qualitatively 
no different from the one used above, provided that (jf x v 
is distinct for each r — 0...Q — 1. However, if v ^ <f) n then 
plaquettes carrying flux v and x v both have the same 
energy cost e m , so that the effective Hamiltonian for the 
spin model now contains dynamical disorder, in the form 
of sites at which there is effectively no transverse field. 

When (f? x v — v for some r, the mapping to Potts 
spins cannot distinguish between pairs of Zq spins x, x+r 
and is hence no longer one-to-one. For example, in the 
SU(2) 2 x SU{2) 2 example of the previous section, pla- 
quettes with flux I are unaffected by the addition of a 
spin-1 vortex. These thus effectively act like (dynamic) 
dilutions in the lattice, meaning that on some sites the 
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ferromagnetic Potts interaction is always satisfied and 
thus adds a constant to the overall energy, independent 
of the spin configuration. (In the more general case these 
dilutions become sites where the NN interaction is satis- 
fied provided that the Zq spins differ by any multiple of 
r) . This is reminiscent of the effect of annealed vacancies 
on the classical 31? Potts model, which has been studied 
in the context of the Blume-Emery-Griffiths^ model. In 
this case, for small perturbations the vacancies are ir- 
relevant at the Ising critical poinlP*', but can drive the 
transition first order at larger vacancy concentrations. 
We therefore expect this to be an irrelevant perturbation 
at the Ising fixed point, and that it will not alter the 
first-order character of the Potts transition. 

When electric sources for the confined edge labels are 
present in the ground state, the mapping to the Zq spin 
model necessarily breaks down. These labels are associ- 



ated with domain walls in the ferromagnetic Potts phase; 
including sources for these labels amounts in the Potts 
language to having open domain walls, which is impos- 
sible in the purely statistical mechanical picture. To in- 
clude these excitations we must include gauge fields in 
the dual theory, which account for the phase winding of 
branch cut singularities at the end of each domain wall. 
In the presence of gauge fields the Landau framework no 
longer strictly applies as the order parameter is necessar- 
ily non-local. 

When H( E ' contains source terms for confined exci- 
tations, therefore, we must explicitly consider whether 
these source terms are relevant at the critical point. To 
do this, we consider the effect of such source terms on 
expectation values of local operators in perturbation the- 
ory. Specifically, we may expand the ground-state wave- 
function in powers of e according to: 



|*> = |* )+E 



Ei — E 



(*<|ir^|*o)|*i) + X; 



(Ei-EoXEj-Eo) 



(42) 



The expectation value of an operator O in the ground state (42 1 can be evaluated to a specified order in e: 



E 



(Ei - Eo^Ej - E ) 



(*o|-ff (B) |^)(^|0|%)(%|-ff (£) |*o) 



+ . 



(43) 



Here we imagine working in a geometry where the 
ground-state is unique; we will return to the more gen- 
eral case presently. The specific form of the higher-order 
terms will not be important for the qualitative arguments 
we present here; the germane point is that the order e n 
term in Eq. (43) contains n powers of the unperturbed 



source i along some curve C in the lattice: 



electric source gap ey in the denominator, and n powers 
of in the numerator sandwiched between various 

excited states. 



Since H^ E > creates open electric strings which do not 
exist in the unperturbed ground state, inner products in- 
volving the unperturbed ground state are non-vanishing 
only when the net effect of these applications of is 
to create some number of sources, move them some dis- 
tance along the lattice, and re-anhiliate them. Thus all 
of the non- vanishing terms in ( 43 ) can be expressed in 
terms of closed Wilson loop operators W% for a confined 



m\om 



E 



E 



l c u ...c r 



<*o|w£ 1 . 



{h,...in} {Cl,...C n } 



(44) 

where oOq' " l " c are coefficients which must be determined 
by the perturbation theory. Here n is determined by 
the order in perturbation theory to which the result will 
be computed. A term where the total length of all the 
Wilson lines is length m can only occur at order m or 
higher in perturbation theory. 

Since the Wilson loop operators Wq can be mapped 
exactly onto operators in the Zq spin model, we may 
study the fate of the critical theory by asking what effect 
these have on the expectation of the spin model's Landau 
free energy. This amounts to asking whether the Wilson 
loop maps to a relevant or irrelevant operator near the 
critical point. 

To answer this question, we must first consider in detail 
the form of the operators in the Zq spin description. We 
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will assume that contains only local terms, which 

either pair-create or move sources within some fixed ra- 
dius on the lattice. This implies that longer Wilson lines 
are suppressed by higher powers of e/ey. A single non 
self-intersecting Wilson line W l c maps to an operator that 
flips a set of Zq spins to generate a domain wall in which 
the Zq spin changes by i across the curve C on the dual 
lattice. Likewise multiple non-intersecting Wilson lines 
map to an operator creating multiple domains. To deal 
with intersecting Wilson lines, we must account not only 
for the positions of the domain walls, but also in general 
for an extra phase that occurs in the gauge theory when 
two Wilson lines cross. This phase is dictated by the 
topological properties of the Wilson lines, and depends 
on the ordering of the Wilson line operators, as illus- 
trated in Fig. [3j If H( E ) creates only a small density of 
defects, then such crossings can occur only at relatively 
high orders in perturbation theory - and are consequently 
suppressed by a high power of e/ey 



In Eq. ( 42 ) we used non-degenerate perturbation the- 
ory, ignoring the topological ground-state degeneracy. 
This is an appropriate starting point in general since the 
ground states in any case do not mix at low orders in per- 
turbation theory. For any finite-sized system, however, 
there will be some order in perturbation theory at which 
it becomes possible to create a Wilson line which wraps 
around one of the non-contractible curves (where these 
exist). For confined sources this maps the Zq spin model 
into a sector with twisted boundary conditions; for de- 
confincd sources it simply introduces a different ground 
state sector into the problem. Here we assume that for 
local operators such effects occur at sufficiently high or- 
der in e/ey that they do not play an important role in 
the physics; indeed if they did, the ground state degen- 
eracy of the topological phase would not be robust to 
their presence, indicating that the perturbation has al- 
tered the underlying topological order before H^ 1 ' tunes 
the system through the critical point of interest to us. 

Restricting ourselves to orders in perturbation theory 
at which such operators cannot occur, we may assess the 
importance of H^ E ' at the critical point by considering 
whether the operator flipping clusters of spins is rele- 
vant. When is local, at finite order in perturbation 
theory the maximum size of the clusters to be flipped is 
finite, and hence the operator is local at sufficiently long 
wavelengths, so that the Landau paradigm applies. In 
the Ising model, where the phase transition is second or- 
der, the Wilson line operators which flip clusters of spins 
are irrelevant at the critical point, and hence generically 
we expect that is as well. This is in agreement 
with the results of Fradkin and ShenkeiHSI, and has been 
verified numerically for the Ising transition between the 
Toric code and the vacuum by Ref. |60j who established 
that electrical sources are irrelevant to the critical theory 
up to the point where the Z 2 topological order itself is 
destroyed. 




(b) 



FIG. 3: Mapping intersecting Wilson line operators to Zq 
domain walls, (a) We specify a set of Wilson lines by their 
contour on the lattice and their linking. Here the red line 
crosses once over and once under the blue line, so that the 
Wilson lines are linked, (b) Up to constants, which are scale 
independent and determined by the topological order, the ef- 
fect of a pair of linked Wilson lines is to flip all Zq spins encir- 
cled by each Wilson line by an appropriate amount. (Here we 
show this for Q = 2, where there is only one type of domain 
wall). Spins encircled by two domain walls will be flipped 
twice. 



Thus the fate of the phase diagram described in the 
previous section in the presence of generic perturbations 
is as follows. The topological order of the gapped phases 
is unchanged so long as the perturbation does not close 
the excitation gap. If the phase transition is second or- 
der, we may evaluate the effect of perturbations on the 
critical theory by leveraging the Landau-Ginzburg theory 
of the corresponding phase transition in the spin model, 
and considering whether the dual perturbation is relevant 
at the critical fixed point. A small density of dynamic 
electric sources, which is dual to flipping clusters of spins, 
gives an irrelevant operator in the Zg spin description. 
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A small density of dynamic magnetic sources is dual to 
annealed disorder in the corresponding 3D classical spin 
model, and is also not expected to alter the nature of the 
phase transition. 



of some excitations into multiple distinct quasi-particle 
types. Here we will describe how these effects arise in the 
condensed phase of the lattice model, explaining how the 
spectrum predicted by Ref. [3T] arises in practice from the 
excitations of the initial Levin- Wen model. 



IV. TOPOLOGICAL ORDER OF THE 
CONDENSED PHASE 



A. String operators and excitations 



In Sect. |TT]we described a lattice Hamiltonian which we 
could tune exactly through a phase transition in which a 
magnetic bosonic excitation condensed. We showed that 
the long-wavelength behavior near the phase transition 
could be mapped onto that of a ferromagnetic Zq spin 
model, and that in the condensed phase certain string la- 
bels - corresponding to domain walls in the spin descrip- 
tion - became confined. In fact, we found that there is a 
special point in the confined phase at which the Hamil- 
tonian Q again becomes exactly solvable, as we may 
consistently project out the edge labels which become 
confined. Here we study in more detail the physics of the 
confined phase, which can be understood by studying this 
second solvable point in the phase diagram. Our objec- 
tives in doing this are two-fold: first, we will see how the 
features of the final-state spectrum (identified by Ref. I3T|) 
arise in the lattice model, and discuss the explicit form of 
the final-state quasi-particle operators. Second, we com- 
ment on the general structure of the lattice models that 
can be obtained as condensates of doubled Chern-Simons 
theories. 

To understand the physics of the condensed phase, it 
is useful to consider the topological properties of an s- 
wave superconductor^. Before the onset of supercon- 
ductivity (in the uncondensed phase), the system is well 
described in terms of electrons and holes (or Fermi liq- 
uid quasi-particles with charge ±e), and the electromag- 
netic gauge field. In the superconducting phase, the low- 
energy degrees of freedom are the Bogoliubov-deGennes 
quasi-particles and the superconducting vortex of flux £ . 
Electrons and holes of the original Fermi liquid theory 
are indistinguishable in the superconductor, since they 
are mixed by the condensate. Further, the Mcissner ef- 
fect confines any gauge field flux unless it is appropriately 
quantized in units of This is necessary so that the 
condensate is single- valued as it winds around the vortex; 
smaller flux quanta would necessarily result in a costly 
branch cut in the condensate wave-function. 

Returning to the question of more general topologi- 
cal symmetry breaking transitions, a general prescrip- 
tion for obtaining the topological properties of the con- 
densed phase is given by Bais and SlingerlancP^, who 
identify three effects of condensation on the excitations 
of the original model. These consist of confinement of 
excitations which braid non-trivially with the conden- 
sate (the Meissner effect), identification of pairs of ex- 
citations which are mixed by their interactions with the 
condensate (analogous to mixing of electrons and holes 
via scattering from Cooper pairs), and possible splitting 
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FIG. 4: Quasi-particle creation operators in the Levin- Wen 
model are 'string' operators which act on states so as to create 
excitations at each of the string's endpoints. (a) These strings 
can carry electric flux, in which case the string operator raises 
or lowers the label on each edge it runs parallel to, as shown. 
In this figure a string operator with label s acts on an edge 
carrying flux i to form an edge with flux i x s. (b) Strings 
can also carry magnetic flux, in which case the operator as- 
signs a label-dependent phase to the wave function each time 
the string crosses an edge separating two plaquettes. In the 
doubled Chern-Simons theory that we begin with, all quasi- 
particle strings are composed of composites of right- (electric 
only) and left- (composite electric and magnetic) components. 
In this figure a magnetic string operator labeled sxs acts on 
an edge i. The net electric flux is thus the combined flux from 
the two s operators, which in general may take on multiple 
values. Thus the vortex string operator is the sum of the 
string operator Ejcresxs- ^or eac h X, Y in the superposi- 
tion the operator acts on the edge by raising its flux by X, 
and incurring a phase, as shown on the right. This operation 
comes with an overall numerical pre-factor not shown here. 
The chief importance of this pre-factor is to ensure that only 
diagrams in which X and Y raise the label i by the same 
amount occur in the superposition. 



Before studying the excitations in the condensed phase, 
we must understand in more detail excitations in the un- 
condensed model. These are most simply described in 
terms of the quasi-particles of the solvable Levin- Wen 
Hamiltonian, which we will describe briefly here. 

As mentioned in the previous section, we may loosely 
speaking divide the excitations of the Hamiltonian ([I]) 
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into constituents which carry electric charge, and those 
which carry magnetic flux. Operators creating electric 
charges on vertices V{ and Vj are open Wilson lines run- 
ning from Vi to Vj (FigureEl). For every edge label i there 
is such a Wilson line, ana consequently a distinct elec- 
tric source. As one would expect from lattice Yang-Mills 
theory, the Wilson line of the source i raises the value of 
the electric flux on the edges that it traverses by i. 

Magnetic sources - such as the $ = 4>r x parti- 
cle described above - are created by an operator that 
essentially assigns a phase to each configuration of edge 
labels. There is one type of magnetic source for each 
allowed edge label i. Unlike in Yang-Mills theory, how- 
ever, in the doubled Chern-Simons theory the elementary 
magnetic source (i^) also carries electric charge^. What 
we refer to as a magnetic i vortex is the achiral particle 
is x ii. When i is a simple current, this is an excitation 
which violates only the plaquette terms of the Hamilto- 
nian, exactly as we expect for a vortex. When i is not a 
simple current, however, its action on an edge e involves 
a superposition of terms, each of which raises e by an 
element of i x i. The precise form of this edge operator 
for a magnetic source is described in detail in Appendix 
[S] and in Ref . HTJ 



B. Confinement and identification in the 
condensed phase 

We begin by understanding how the first two phenom- 
ena - identification of some excitations, and confinement 
of others - occur in the condensed phase of Q. Es- 
sentially, this is no different than in the superconductor. 
First, as we have already seen, excitations which braid 
non-trivially with the condensate (those for which qi ^ 0) 
are confined in the condensed phase. This is analogous 
to the Meissner effect, since their braiding statistics with 
the condensate ensure that they engender branch cuts 
(or domain walls, in the Zq spin picture) in the conden- 
sate wave function. Further, after condensation, pairs 
of particles which are mixed by scattering with the con- 
densate are indistinguishable (or identified) in the final 
topological theory, much like electrons and holes in a su- 
perconductor. 

We have already demonstrated confinement explicitly 
in the previous section. Specifically, the arguments used 
to show that closed domain walls incur a linear energy 
cost apply equally well to static sources: though their 
energy can exceed the cost 2ey of spontaneously creat- 
ing a new pair of sources, the number of these sources 
is conserved by the Hamiltonian ^ and such processes 
cannot occur. At the solvable point where edge labels 
with q, L 7^ do not appear in the ground state, the ener- 
getic cost of adding a pair of sources with ^ is linear 
in the length of the string separating them, as there must 
be a string of edges with confined labels which connects 
the pair. Away from the solvable point there is in general 
a small admixture of closed strings of the confined Qi ^ 



labels in the ground state, but as this is relatively small 
(domain walls are confined) the energetic cost of a pair of 
confined sources will still scale linearly with their separa- 
tion. (As in QCD, perturbations mixing confined sources 
(or open confined strings) with the ground state lead to 
difficulties in identifying the confinement potential, as 
the string joining a pair of test sources can be broken 
into shorter segments by creating confined particle-anti- 
particle pairs). 

To show in detail how identification occurs, we note 
that in the ground state of Eq. [7] at ar = 0, ajv = 1 
each plaquette is in a superposition of states with exci- 
tations Id, $, ..<i> ( 5 _1 . Hence any string operator W^ X s 
applied to this phase will result in the same superposition 
of sources with the labels s,$xs, ..$ < 2~ 1 x s. (Here s is a 
particle of the doubled theory). In the condensed phase, 
therefore, there is a single string operator associated with 
the sources 

s, $ x s, ...<i> <3 ~ 1 x s (45) 

and correspondingly all of these should be identified as 
creation operators for the same quasi-particle. 

The rules for mixing here are inherited from the 'fusion' 
rules, which specify how to combine representation labels 
on the edges of the lattice. In a topological theory there 
are rather stringent consistency constraints on these rules 
(see, for example, Ref. 135])) but for the present discussion 
the important point is that the excitations naturally mix 
with the condensate in a specific way. When $ is a simple 
current, the fusion rules specify that $ x s is a single 
quasi-particle species. 

We illustrate the effect of confinement and identifica- 
tion on the spectrum with a few examples. It is conve- 
nient to define the notational convention 

(r, s) = r R x s L (46) 

for excitations in the doubled Chern-simons theory. In 
this notation, the magnetic $ = (f> R x tfi L particle that 
we condense has the form (0, <fi). 

1. SU(2) k fork odd 

The labels appropriate to SU(2)k Chern-Simons the- 
ory are the total spin of representations s — 0(= 
Id), 1/2, ...k/2. We condense the (k/2,k/2) excitation, 
which is an achiral simple current of order 2 (i.e., a Z2 
boson). 

To find the confined particles, we note that the S- 
matrix elements are -^-St s = (— l) 2s (see discussion 
near Eq. [51). Hence edge labels with half-integer spin 
are confined (they all are effectively mapped to domain 
walls of the Ising model), while those with integer spin 
remain deconfined. In the solvable limit, this implies that 
the plaquette projector also contains only integer spin la- 
bels. Since the edge label (or electric flux) depends on the 
combined electric fluxes from the right- and left-handed 
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components of each excitation, the deconfmed particle 
types are: 



(r, s) and 



(£•*) 



(47) 



for r and s integer. 

Fusing the (k/2, k/2) particle with other excitations in 
the theory gives: 



(k/2, k/2) x (r, s) = (k/2 - r, k/2 - s) 



(48) 



As promised, the rules for adding angular momenta in 
the Chern-Simons theory have been deformed such that 
the product on the left-hand side gives a unique result on 
the right. Eq. (48) therefore identifies pairs of excitations 



(r, s) 



2r k-2s 



(49) 



When k is odd, if r and s are integers, then both labels 
on the right are half-integer. Thus we may eliminate 
the half-integer labels from the theory completely, and 
identify each excitation by the appropriate integer labels. 

In this case, this is the whole story: condensation has 
merely eliminated all half-integer spin labels from the 
theory. The result is a doubled SO(3)k Chern-Simons 
theory!^. 



2. SU(2) k for k even 

After condensing the (k/2, k/2) boson, if k is even, our 
result will differ from that described above in that k/2 is 
an integer, so that fusion with the condensate identifies 
integer pairs with integer pairs, and half-integer pairs 
with half-integer pairs. Hence we can no longer eliminate 
the half-integer labels entirely from the theory, as they 
may still appear in the guise of excitations of the form 
(r/2,s/2) with r and s odd, even though they will no 
longer enter as deconfmed edge labels. This suggests that 
the topological description of the condensed phase is not 
a doubled Chern-Simons theory as we shall see further 
below. 

A second peculiarity of the identifications is that the 
quasi-particle (fc/4, k/A) maps to itself under fusion with 
the condensed (k/2, k/2). We will see that this is not 
unrelated to the fact that the quasi-particle spectrum 
cannot be simply separated into right- and left- handed 
components when we discuss splitting. 



C. Splitting 

The final possible effect of condensation is that some 
excitations will split into multiple distinct quasi-particle 
species. (This phenomenon occurs in spontaneously bro- 
ken non-Abelian gauge theories, but not in the super- 
conducting example discussed above). It is important to 



note that this splitting does not change the total dimen- 
sion of the low-energy Hilbert space. Instead, it splits 
a 2- (or more generally d ) dimensional internal Hilbert 
space of the excitation in the uncondensed phase into 
multiple 1- (or d/n) -dimensional Hilbert spaces. 

An intuitive understanding of this splitting can be 
gained by considering how it arises in non-Abelian gauge 
theories. Essentially what happens is that the excitation 
in question transforms under a d-dimensional represen- 
tation of an internal symmetry which is broken in the 
condensed phase into multiple lower-dimensional repre- 
sentations of the residual symmetry group. For exam- 
ple, in a theory with unbroken SU(2) symmetry, a spin- 
1/2 particle is associated with a 2-dimensional internal 
Hilbert space. (In other words, we may express particles 
transforming in the fundamental representation of SU(2) 
as 2-component vectors). If we break the SU(2) symme- 
try by condensing a spin-1- Higgs field, this 2-component 
vector can be separated into its spin-up and spin-down 
constituents, which are now no longer related by sym- 
metry. The residual U(l) gauge transformation acts on 
these as: 



e c 



t J 



(50) 



In other words, the 2-dimensional representation of 
SU(2) has separated into two 1-dimensional representa- 
tions of U(l) (here carrying opposite charges). The total 
dimension of the Hilbert space associated with each par- 
ticle creation operator is still 2, though we now have 
distinguishable spin- up and spin-down excitations. 

One indicator that such a splitting had to happen in 
the above example can be found in the rules for combin- 
ing the representations of SU(2). Combining two spin 
1/2 excitations gives: 



2 2 



1 



(51) 



which is a superposition of the (gauge-neutral) singlet 
and the spin-1 triplet excitations. Before condensation, 
these are distinct excitations, as they transform in differ- 
ent representations of the symmetry group. After con- 
densation, however, the residual U(l) symmetry group 
cannot distinguish between the singlet and triplet states. 
Instead, it is sensitive only to the l z eigenvalue of each 
state, rather than to I. Labelling states on the right-hand 
side of (51) by their l z eigenvalues, we have: 



\ x ^=0 + + l + (-l) 



(52) 



which is to say, on the right-hand side we obtain two 
distinct copies of the singlet (l z — 0) representation of 

,t J 



U(l). (In this case, these are cLcL ± c 



The rules 

of representation theory dictate that two copies of the 
singlet can be obtained on the right-hand side only if 
there are two distinct excitations on the left (which we 
may choose to be cLct). Hence we conclude, solely by 
examining the rules for addition of angular momenta, 



16 



that the spin- 1/2 excitation had to split into two dis- 
tinct l-dimensional representations (and hence 2 distinct 
quasi-particle types) in the condensed phase. 

This same logic about combining representations in the 
condensed phase applies to the more general framework 
of topological symmetry breaking^SHSl] This will indi- 
cate that, in some of the examples of the type discussed 
here, certain excitations split into multiple distinct quasi- 
particles after condensation. We will first review the cri- 
teria for splitting to occur, and then explore how it arises 
in practice in the lattice models. 



1. Determining whether quasi-particles split 

The generalization of the criteria we found above in 
SU(2) Yang-Mills theory for splitting after condensation 
is as follows. Suppose <& is a simple current (with $® = 1) 
which we will condense. Then an excitation r will be split 
in the condensed phase if 



(53) 



for k < Q. As in the example above, the reason for the 
splitting is that in the condensed phase, r x r contains 
multiple copies of the trivial representation, and hence 
must be split into multiple particle types if the represen- 
tation theory is to remain consistent. Specifically, the 
representation f is by definition the one which combines 
with r to give the the singlet representation (plus some 
other representations, in general). Thus we have: 



Id 



r x (r x $ fe ) = $ fe - 
r x (r x $ Q ) = $ Q 



which implies that 
r x r = Id + <& & 



2/,- 



<I><2 



(54) 



(55) 



That is, the tensor product of representations r and r 
contains (among other things) all powers of $ fc (mod Q). 
Since in the condensed phase any power of $ is identified 
with the trivial representation (since, as before, the con- 
densate is by definition in the singlet representation of 
the residual symmetry group) , this gives ? singlet repre- 
sentations on the right-hand side of (55). Excitations in 
representation r before condensation consequently split 
into ^ distinguishable particle types in the condensed 
phase. 

It is a feature of the representation theory that the 
total dimension of the Hilbert space (total quantum di- 
mension) is preserved by this splitting. In particular, if 
r is l-dimensional then it follows that k = Q and the ex- 
citation cannot split. This is rather obvious in the case 
that r is truly a representation of a non-Abelian symme- 
try group; however, it also holds true in the truncated 
representation theory germane to the Chern-Simons lat- 
tice models considered here^. 



2. Splitting on the lattice 

Armed with this simple criterion to understand when 
some excitations in the condensed phase will split, we 
now turn to the question of how this splitting manifests 
itself on the lattice. In the example given above, we 
could explicitly identify the two split spin- 1/2 particles 
as we knew precisely the form of the residual symmetry 
generator's action on the members of the original SU(2) 
multiplet. In the lattice model Q we do not have ac- 
cess to this information. However, we will be able to 
identify a set of labels (analogous to S z in the example 
above) which are indeterminate before condensation but 
separately conserved in the condensed phase. 

Though this may seem like a rather trivial exercise in 
practice, since we have already argued on general grounds 
that such splittings must occur, it is actually important 
to demonstrate that splitting occurs in order to conclude 
that the condensed phase does indeed represent a con- 
sistent topological theory. As we discuss briefly in the 
conclusions, there are situations where it is not clear that 
this is the case. 

Let us begin with an example, and consider condensing 
the spin-1 excitation in doubled SU(2) 2 Chern-Simons 
theory. The fusion rules for this theory are given in 



Eq. ( 26 ) , and in particular stipulate that 



(1,1) x 











a 




■a 





(56) 



so that the achiral spin- 1/2 particle must split into 2 
distinct excitations after condensation of the (1,1) bo- 
son. To identify these distinct excitations, we first note 
that the achiral spin- 1/2 particle carries a magnetic flux 
from its 1/2/, component, and an electric flux from the 
combination of its 1/2^ and 1/2^ components. Since 
1/2 x 1/2 = + 1, the electric flux associated with this 
excitation can be either or 1 on a particular edge. The 
precise form of the string operator dictates that it may 
change between and 1 when the string crosses between 
two plaquettes over an edge carrying a spin- 1/2 label 
(see Appendix [A]). At the solvable point in the condensed 
phase, however, the spin-1/2 edge labels have been com- 
pletely eliminated from the Hilbert space. Thus in this 
limit, a (|, |) particle is associated with an electric flux 
which is either or 1 at all points along the string op- 
erator (and consequently, also at the vertices on which 
the string terminates). We may therefore identify two 
distinct quasi-particle types, (|, |) and (|, 

The key point here is that prior to condensation, a 
string which is purely of the or (5, 5), type is 

not topological. That is, even at the exactly solvable 
point this string is physically observable - whereas the 
string operators of Ref. 36] create strings for which only 
the end-points have physical meaning. The (|, |) Q and 

(2' \)i strings, however, create a spin-1/2 vortex at each 
endpoint, together with a string of plaquettes in a super- 
position of the ground state and the spin 1-vortex (the 
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(1, 1) boson) excited along its trajectory, as explained in 
Appendix lAl Hence in the uncondensed phase, the eigen- 
statc ( | , |j is associated with a 2-dimensional Hilbcrt 
space (of electric flux or 1). In the condensed phase, 
where the string of possible spin 1-vortex excitations is 
undetectable, (|, |) and (|, |) 1 are two distinct topo- 
logical quasi-particle creation operators. 

This example, though relatively simple, illustrates pre- 
cisely how splitting occurs in general condensates. First, 
notice that if (</> fe , <p k ) x (r, s) — (r, s), then we must have 
both 4> k x r — r and <fi x s = s. The electric flux asso- 
ciated with this excitation is a superposition of 

rx s= I . (57) 

igrxs 

After condensation, we find that the sum on the right- 
hand side will split into subsets of labels which do not 
mix once the confined edge labels have been projected 
out of the theory. It is possible to deduce from the fact 
that $ fc x (r,s) = (r, s), and that $ acts as the iden- 
tity in the condensed phase, that there are § excitations 
which are topological in the condensed phase. (As above, 
in the uncondensed phase there is only one). A (rather 
technical) proof of this fact is given in Appendix [Bj 

Thus we find that the excitations in the condensed 
phase split into § distinct species of excitation, as re- 
quired for a consistent topological phase. 



left-handed labels. Specifically, any deconfined excitation 
(a, b) has q a = lb- Since charge is additive under fusion 
with cj) k , and since by assumption the charge of 4> k spans 
all possible values, there is some r for which q$r = —q a . 
Hence (<ff ,<t> r ) X (a,b) is composed of two string-types 
which both have g-charge 0. Hence every excitation in 
these theories can be viewed as a composite of two decon- 
fined particle types. In this case, the condensed theory 
is just two opposite chirality copies of a subset of the 
particles of the original chiral Chern-Simons theory. 

If q^k = for some k < Q, however, the above result 
need not hold. It is easy to construct examples of this 
in Abelian Chern-Simons theories. For example, the k 
particle in U(l)2k has: 

kxk = Id S kj = (-l) j . (58) 

If k is odd, then all deconfined excitations can be ex- 
pressed in terms of pairs (2j, 21) of deconfined excitations, 
and the theory is again a tensor product U(l)kX U(l) k of 
identical right- and left- handed Abelian Chern-Simons 
theories. If k is even, the deconfined excitations fall into 
two classes: (2j, 21) and (2j + 1,21 + 1) which are not 
equivalent under fusion with (k, k). (When k is odd these 
odd and even sectors are identified) . Hence here the spec- 
trum is not a direct product of two chiral components, 
as neither component of the odd excitations can exist in 
isolation. 



D. Structure of the final theory: Examples 



2. SU(2) 2 



Having established the nature of the spectrum in the 
condensed phase, it is worth pausing to take stock of the 
variety of possible topological phases which can be cre- 
ated in this way by condensing a magnetic simple current 
in a doubled Chern-Simons lattice model. We will illus- 
trate this with a series of examples. 



Theories without splitting 



The simplest case where splitting does occur is after 
condensation of the achiral spin 1 excitation in a dou- 
bled SU {2)2 Chern-Simons theory^!. As noted above, 
the (1/2,1/2) excitation splits into two components in 
the condensed phase, distinguished by their electric flux 
(which may be that of the spin singlet or spin triplet). 
Since the chiral spin-1/2 excitations are confined, this 
leaves us with the following 3 quasi-particles in the con- 
densed phase: 



The simplest case we may consider is that of a the- 
ory in which no particles split in the condensed phase. 
The properties of the spectrum here are determined by 
confinement and identification alone. Depending on the 
nature of the condensed excitation, the final theory may 
be a new doubled Chern-Simons theory (with a gauge 
group that is a quotient group of the original, as gen- 
erally occurs when vortices are condensed^) . Alterna- 
tively, it may be a theory in which not all excitations 
can be decomposed into separate right- and left- handed 
chiral components. 

The difference between these two cases is determined 

by S^k^ = e Q , where we condense (</>, <j)). If q^k ^ 
for all k < Q, then in the condensed phase each set 
of identified quasi-particles either is confined or contains 
one element which is a composite of deconfined right- and 



(0,1) = (1,0) 













2)0 







(59) 



The first of these, which is a purely electric source in our 
lattice model, is a fermion. The other two (which we 
identify as purely magnetic, and both magnetic and elec- 
tric) are bosons (but here with relative semionic statis- 
tics). 

These excitations give precisely the spectrum of Ki- 
taev's Toric code^ (or Z2 gauge theory with matter) . In- 
deed, at the solvable point in the condensed phase, where 
we eliminate all spin-1/2 edge labels from the theory, we 
may use the edge-labeling scheme: 



(60) 
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In this basis, and dropping the terms —(—1)™" (which we 
take always to be —1, since n a = 0) the Hamiltonian is 
precisely that of the Toric code: 

tf = -£IK-£II-e (61) 

V e P e 

In matching the excitation spectra to that of the Toric 
code, we must contend with one subtlety of the condensed 
phase - namely, the purely electric source is fermionic, 
rather than bosonic as it should be for the Zi gauge the- 
ory. The reason is that operator which creates the (1, 0) 
excitation is an electric-type string which raises the spin 
on each edge by 1 (mod 2); however, it also obtains a 
phase of for each edge e it crosses. The operator cre- 
ating j). excitations is a magnetic-type string which 
assigns a phase of —1 for each edge of spin 1 (mod 2). 
(I' l)i * s ^ ne operator which raises the spin by 1 without 
inducing any phases - and is mutually semionic relative 
to both (1,0) and (|, |) . In the spin basis, this gives 
the quasi-particlc operators the final form: 

•s(i/2,i/2)o = ny* %/2,i/2)i = ny* 

e e 
S(1,0) = 5(1/2,1/2)05(1/2,1/2)! ( 62 ) 

which identifies (|, |) as the electric source of the Zi 
gauge theory, (~, ~) as the Z2 magnetic source, and 
(1, 0) as their (fermionic) composite. 

3. SU(2) k for k even 

We may generalize some of the features of the SU(2)2 
example above to SU{2)k for general even k, condensing 
the (k/2,k/2) vortex. The deconfined excitations here 
have net integer spin on each link, and hence must have 
the form with i and j either both integer or both 

half-integer spins. This gives 2 (|) + 2| + 1 excitations 
before identification. Of these excitations, all but (f , f ) 
get identified in pairs; (f, 4) in fact splits into two excita- 
tions. Thus we obtain a total of (|) 2 + | + 2 excitations. 
For k > 2 the number of excitations is generally not a 
perfect square, so that the topological order cannot be 
that of a doubled Chern-Simons theory (or any double 
T x T of a valid topological theory) . 

When k = 2, we showed above that the condensed 
phase is the topological limit of a discrete gauge theory. 
For k — 4, the condensed phase is also described by a 
discrete gauge theory, in this case a twisted version of 
the non-abelian gauge group D3. For k > 4, however, the 
topological order of the condensed phase cannot be that 
of a discrete gauge theory, since some of the deconfined 
particles have non-integer quantum dimensions. 

We can nonetheless write the complete set of string 
operators for these theories. There is an even sector of 
string operators which are composed entirely of the inte- 
ger spins of the uncondensed phase, and an odd sector of 



string operators composed of pairs il* Jr of half- integer 
spins of the uncondensed theory. The even sector con- 
sists of excitations which can be constructed using only 
strings which correspond to deconfined labels. Specifi- 
cally, we have: 

. , \iL,in (i = /<£..§- 1) 

simple : < k f 
I I (non-chiral) 

composite : ir x Jl i < f , j < § (63) 

with i and j integer. Because of the identifications (which 
identify even excitations with even, and odd with odd, 
unlike the scenario for odd k), the excitations in the 
even sector contain only some of the allowed composites 
which we can construct from the chiral string operators 
%l,r- In addition, there is only a single string operator 
\ L = \ R associated with the | particle, whose chirality 
is no longer well-defined after condensation. Indeed, any 
particle of the form xr x ( | — x) L is non-chiral in the 
sense that in the condensed phase it is indistinguishable 
from the opposite chirality excitation xt, x (| — x ) r/ 

In addition to these, there are excitations which can- 
not be composed of simple string operators. That is, they 
are composites of pairs of half-integer spin string oper- 
ators. These composites have become irreducible in the 
condensed phase, where the individual half-integer spin 
strings have been confined. These excitations are: 

k k 

iR*jL i<^,j<2 (64) 

with i and j half-integer. 

The split (2, 2) particle is in the even sector for k = 
(mod 4), and the odd sector for k — 2 (mod 4). 

4- Drinfeld Doubling 

The general structure of the construction here is as 
follows. At the solvable point in the condensed phase, 
we may project onto states composed only of deconfined 
edge labels (those which braid trivially with the con- 
densed magnetic excitation), from which our fixed-point 
Hamiltonian for the condensed phase (Eq. [m] ) and its 
low-energy excited states are constructed. When q^k 7^ 
for any k < Q, all excitations in the final theory can be 
constructed from string operators containing only these 
deconfined edge labels. When q^k = for some k < Q 
(which must occur if there is splitting, but may occur 
in other examples as well, such as the Abelian theories 
discussed above), we generally find that not all excita- 
tions in the final model can be expressed in terms of 
strings corresponding to deconfined edge labels. Never- 
theless, the resulting topological theory is equivalent to 
a Levin- Wen model built on the category of deconfined 
edge labels only. In general a Levin- Wen model built 
from a category produces a topological theory known as 
the Drinfeld double of the category. In our case we iden- 
tify the "non-simple" string operators of Ref. |3S]as those 
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which cannot be constructed from deconfined edge labels 
alone, but require additional phases to account for the 
fact that they are composites of pairs of confined electric 
sources. 

The simple protocol outlined here shows how certain 
Drinfcld doubles can be interpreted physically as the out- 
come of condensation in a Chern-Simons theory. Many of 
the examples given above are in any case discrete gauge 
theories, which we could alternatively understand as the 
result of Higgsing a continuous Maxwell or Yang-Mills 
theory. However, some (for example, those obtained 
by condensing the (k,k) particle in a doubled S , t/(2) 2 fc 
Chern-Simons theory with k > 2) are not; this approach 
offers a physical mechanism for the origin of these states 
from models which can at least be understood in terms of 
continuous field theories. In general, our approach sug- 
gests that we may view some Drinfeld doubles as aris- 
ing because an achiral excitation condenses in a doubled 
Chern-Simons theory. Specifically, because the conden- 
sate is achiral, time reversal symmetry is preserved on 
both sides of the phase transition, though the spectrum 
of the final theory cannot always be decomposed into de- 
coupled right- and left- chiral sectors. 



V. CONCLUSIONS 

In this work we have given an explicit realization of 
topological symmetry breaking^ in lattice models. By 
constructing a lattice Hamiltonian that can be tuned be- 
tween two solvable Levin- Wen points, via a condensation 
transition, we may map the topological symmetry break- 
ing transition explicitly onto a 2 + ID transverse-field 
Potts transition. The phase transition can be under- 
stood by studying the dual Potts description, to which 
the Landau formalism applies if the transition is second 
order. Though this duality is precisely valid only for a 
very special trajectory through the phase diagram, we 
argue perturbatively that the effect of small deviations 
from this trajectory can be understood within the Lan- 
dau theory of the spin model, and will be irrelevant at the 
critical point. This gives a general framework to clarify 
the relationship between phase transitions separating dif- 
ferent topological orders (such as those described by ref. 
I34j) and phase transitions of the Landau type. Further, 
we have studied the properties of the condensed phase, 
and identify the complete set of quasi-particle creation 
operators required for a consistent topological phase. 

The type of transitions we discuss here are special in 
two ways. First, in topological theories, two bosons of 
the same type generally combine to give a variety of other 
species of bosons. This is analogous to combining spins, 
where for example 1/2 x 1/2 = 001, and occurs be- 
cause, in the same way that particles can be classified 
by their transformation properties under rotations (or 
total spin) , excitations in a topological theory are associ- 
ated with representations of a (quantum) grourJ^. When 
the group in question is non-abelian, most excitations 



will not be simple currents. The general technique em- 
ployed here to construct the Hamiltonian Q by adding 
a term that pair-creates vortices on adjacent plaquettes 
still applies in such cases. Condensing bosons with non- 
Abelian fusion rules will lead to a different critical be- 
havior, which is not equivalent to that of any statistical 
mechanical model that the authors are aware of. The 
study of these transitions is undoubtedly a rich subject 
for further study; one interesting example is discussed in 
Ref. [H 

The second restriction we have imposed here is to con- 
sider only achiral condensates, by condensing plaquette 
violations in the lattice model. An obvious question is 
whether chiral condensates (or condensates of vertex vi- 
olations, in the lattice model) can also occur. From the 
purely topological viewpoint there is no obstruction to 
forming these^S, provided that the excitation to be con- 
densed is a boson. (This is always the case for achiral 
excitations, but need not be for their chiral cousins) . Fur- 
ther, the critical theory will again be of the transverse- 
field Potts type if the condensed boson is a simple cur- 
rent. In the lattice model, however, only when there is 
no splitting is it clear that operators for all excitations 
in the condensed phase can be constructed. In the ab- 
sence of splitting, the Chain-MaiP^ formulation of the 
partition function can be used to show that the final the- 
ory is dual to an achiral condensate, and the resulting 
duality mapping between the string operators gives an 
explicit representation of all excitations in the final the- 
ory When splitting occurs, this duality fails and there 
appear to be no conserved quantum numbers to differ- 
entiate the split particle types, suggesting that the final 
topological phase may not be fully realizeable by the lat- 
tice model. We will discuss these results in more detail 
in a future work. 

The task of fully categorizing the possible phase tran- 
sitions and critical theories between phases of different 
topological order remains a source of many open ques- 
tions. The solvable Levin- WerP^ models considered here 
provide a useful framework in which to rigorously study 
such questions; since both topological order and prop- 
erties of the critical theory are relatively universal, con- 
clusions drawn from the lattice model also apply to real 
physical systems exhibiting the desired topological char- 
acteristics, where these exist. 



Appendix A: Magnetic quasi-particle operators 

Here we discuss in more detail the form of operators 
creating magnetic quasi-particles which do not corre- 
spond to simple currents. The precise form of the general 
magnetic quasi-particle operators is complicated by the 
fact that, in the uncondensed phase, if the label i does 
not correspond to a simple current, there is no purely 
magnetic excitation associated with i. Rather, what we 
will call the magnetic i excitation is in fact a specific su- 
perposition of excitations which all carry magnetic flux 



i, but also carry electric flux j G i x i. 

The reason for this is that the fundamental excitations 
in our Levin- Wen model are not electric and magnetic 
fluxes as would be the case for a lattice gauge theory, but 
the sources of the right- and left- handed Chern-Simons 
fields. As explained in Ref. STJ the right-handed sources 
are precisely the electric sources described above. The 
left- handed sources i^, however, carry both the electric 
charge i and the magnetic flux i. The best approxima- 
tion to a purely magnetic excitation in this case is the 
achiral source X in, which has magnetic flux i and 
electric flux i X i = Id + .... In general the individual 
electric flux labels on the right-hand side are not con- 
served along the length of a particular string, so that the 
entire superposition is required to construct the appro- 
priate quasi-particle operator. 
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(Id,ip) 


(ip,Id) 
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V> 
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TABLE I: Action of the four possible combinations of fusion 
channels for the ox X or excitation (the a vortex). Here we 
have omitted any factors associated with labels on adjacent 
edges. 

of the link respectively. This action can be decomposed 
into the action of the four channels shown in the figure, 
as shown in Table llj The coefficients Fl^r come from the 
action of the ip tail on the left- and right- vertices. 





FIG. 5: The a vortex in the doubled Ising anyon theory 
consists of a pair of electric sources (one right- and one left- 
handed). Its action on an edge is given by a phase (depicted 
here by the red ring) every time the string crosses between pla- 
quettes, together with an electric component which raises the 
edge label by Id or ip. This gives the four possibilities shown 
here each time the string operator crosses between plaque- 
ttes. The labels (Id, Id), (if), Id), (Id, ip), and (ip,ip) denote 
the associated electric flux on the upper and lower sides of 
the crossing. 



To illustrate how this works in practice, we consider 
the a vortex in the doubled Ising anyon theory. (This is 
essentially the same as the spin-1/2 vortex in a doubled 
SU(2) 2 Chern-Simons theory, though the signs in this 
case are somewhat simpler to track). The operator is 
represented diagramatically in Fig. [5] it acts on the edge 
labeled i in the Figure according to 

3(1/3,1/3)1*"} = \ {Sl /2 ,i\i) + (-l) n ^ 1/2>i |i X V>» 

(Al) 

where F^r = ±1 are coefficients which depend on the 
labels of adjacent edges at the left (L) and right (R) ends 



It is important to note that if we keep only the (Id, Id) 
and/or (ip,ip) fusion channels, the action of this opera- 
tor on a given edge is unchanged if we simultaneously 
act with the operator (—1)™" which creates a pair of ip 
vortices. This means that keeping only the (Id, Id) and 
(ip,ip) channels adds an indefinite number of ip vortices 
to each plaquette. In the condensed phase, where ip vor- 
tices are in any case not conserved, this does not affect 
the energy of states this operator creates; in the uncon- 
densed phase, however, if the a vortices at the end-points 
of the string are more than one plaquette apart then the 
state is not a single eigenstate, but rather a superposition 
over all eigenstates with some number of ip vortices on 
the intervening plaquettes. As a consequence this opera- 
tor is not topological, since the location of these possible 
ip vortices on the lattice depends on its trajectory. Once 
we include the effect of the (Id,ip) and (ip,Id) channels, 
which do not annihilate a labels, this is no longer the 
case and the operator creates only a pair of a vortices at 
each of its endpoints. 

The case for general magnetic sources is similar: we 
find that only when all of the fusion channels of i X i 
are included as electric source lines is the string operator 
topological. If some of these are omitted, then the op- 
erator has some probability of creating extra vortices on 
the plaquettes separating the two i vortices, and hence is 
not topological as the number of possible violations scales 
with the separation between the string's endpoints. 



Appendix B: Splitting in general condensates 

Having detailed the form of composite operators ajx 
an in the previous section, we now present a proof that 
condensing an achiral simple current will always lead to 
the correct splitting of string operators in the condensed 
phase. It is useful to begin with the example discussed 
above, and consider the splitting of the er^ x particle 
in the doubled Ising theory. 

From the action of the possible fusion channels of the 
a vortex in Table [TJ it is easy to see that when the edge 
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label a is eliminated in the topological limit of the con- 
densed phase, the (Id, Id) and (tp,^) channels do not 
mix. One way to understand why this happens is to 
note that before condensation these string operators were 
non-topological because they created a trail of plaquettes 
containing a superposition of no vortex and the "0 vortex; 
hence their average energy depends linearly on the string 
length in this regime. After condensation, however, the 
presence of a ip vortex does not change the energy of 
a state, and these become legitimate topological string 
operators. 

This basic argument can be generalized to other con- 
densates of achiral simple currents. To do so, we ex- 
ploit the fact that topological excitations of the (solv- 
able) Levin- Wen Hamiltonian describing the condensed 
phase must obey the hexagon relations. From the point 
of view of the lattice model, these simply state that any 
string operator creating a topological excitation can slide 
freely over vertices, as shown in Fig. [6] As in the case 
of the olo~r particle, string operators which do not sat- 
isfy the hexagon relation are not topological because they 
leave a trail of excitations along their length, rather than 
just at their endpoints. In some cases, however, this trail 
of excitations consists entirely of the vortices that we 
condense. When this happens the corresponding string 
operators are topological in the condensed phase. Our 
purpose here will be to show that this process accounts 
exactly for the splitting expected from the TSB criterion. 
Readers should note that to do this we will make refer- 



ence to tensors F (the 6j- symbols) and R (the universal 
R matrix) which are determined by the choice of topolog- 
ical order. We will not explain their meaning here, but a 
useful introduction can be found in Refs. I55H41I62"! 




k k 



FIG. 6: String operators creating topological excitations can 
slide freely over vertices, without altering the state of the sys- 
tem. String operators with this property are said to obey the 
hexagon relatioi%^\ the corresponding excitations are topo- 
logical because their energy is independent of the trajectory 
which the string takes through the lattice. String operators 
that are not topological leave behind a trail of excitations. 
After condensation, some achiral particles may split - mean- 
ing that sub-sets of the possible edge labels associated with a 
particle aL X an separately obey the hexagon relation. This 
happens because the trail of excitations that these string op- 
erators leave behind are vortices which have condensed. 



Before condensation, in the doubled Chern-Simons theory all string operators have the form d L x Sr. These can 
be expressed in the form of a phase operator R a b acting on each edge that the string crosses (diagramatically a ring 
labeled (a, b) encircling the crossed edges, as shown in Fig. [5]), and operators sx$y which raise or lower the electric 
flux on each edge by bi and bf to either side of the ring. We will focus on the configuration shown in Fig. [6j in which 
the string operator crosses over a pair of edges. For example, d L x clr has the form: 

d R ®d L =R aa ^Ff,$b . (Bl) 

where F£,£® * s a coefficient (6j symbol) dependent on the labels a and b. The reason that this particle satisfies the 
hexagon relation is that since the sum runs over all values of b allowed by fusion, we may 'undo' the process of making 
the composite operator, to depict this operator as two separate strings, labeled ur and d L respectively. In the doubled 
Chern-Simons models, it is easy to show that the string operators ul, or do obey the hexagon relation^. Graphically, 
we may depict the situation as follows: 




Here X indicates that each a loop is encircling a labeled edge in the lattice, as is appropriate for the phase operator 
R aa . After re-expressing the operator locally in terms of the two strings d L and a#, the resulting operator may be 
pulled over a vertex at which the two edges X join. This ensures that the hexagon relation (Fig. [6]) is obeyed. 

Before condensation, only when the dashed line carries the label can we pull the a-loop on the right-hand side 
over a vertex, and hence there is only one topological string operator associated with d L x or. After condensation, 
however, any combination of coefficients on the left-hand side which results in the dashed line carrying powers of (f> 
gives an operator that satisfies the hexagon relation, and hence a valid topological quasi-particle. Here we assume 
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that there are no other labels in the category which braid trivially with all of the deconfined edge labels; hence we 
wish to find all linearly independent sets of coefficients C(b) such that: 




where the coefficients a n are arbitrary. 



Our task now is to count the number of possible lin- 
early independent sets of coefficients C(b) will have the 
property that 



E 



C(b)F ( 



aa^b* 
aa* d 



(B2) 



n=0 



This will give us the number of different particle types, as 
each dimension in the vector space of possible solutions to 
(B2 1. Clearly, if <p n x a ^ a, then a n — as the diagram 



on the right is not consistent with the fusion rules of the 
theory. This leaves Q/k possible linear combinations on 
the right-hand side of Eq. ( |B2[ ), corresponding to the Q/k 
independent particle types. 

To show this explicitly, we will identify the Q/k lin- 
early independent sets of coefficients C(6). We begin with 
the choice C(b) = ^ = F£* Q a 6 °. With this choice 
have: 



we 



A 





a a 





E b C k (b) 



with b 



C k {b) = 





a'*bb* 



R 



<j> k b 



(B3) 



We thus have a candidate choice of C nk for each of the 
Q/k possible values of n. These must all be linearly in- 
dependent, since the vertical lines of the diagrams on 
the right carry different powers of <fi. Any other choice 
of coefficients either gives a superposition of these Q/k 
possibilities, vanishes, or produces a quasi-particle type 
that does not obey the hexagon relation. In particular, 
a x a contains (jp only if <$P x a = a, so that no other 
powers of <p may appear on the right-hand side. 

It is useful to express this statement in matrix form. 
We can write 



Co (61) 



C (b 2 ) 
C\(b 2 ) 



C (b r ) 

Cx(b r ) 



and hence C (b) = is certainly one of the allowed sets 
of coefficients. Next, we consider diagrams of the form: 





If <j) k x a = a, then the four external edges of the diagram 
are still labeled a. Further, the vertical line on the right- 
hand side clearly carries the label (f> k , and hence the parti- 
cle depicted obeys the hexagon relation in the condensed 
phase. Specifically, the identity above is equivalent to: 



Cq/k{b\) C Q / k {b 2 ) ... C Q / k (b r ) 

where r is the number of possible fusion outcomes of a£g)a, 
and in particular r > Q/k as a ® a = Id + <j> > + ... + 
4> N ~ k + .... The Q/k orthogonal linear combinations of 
coefficient vectors which can be formed from these ensure 
the existence of Q/k distinct string operators. In general 
these string operators will generate linear combinations 
of the true quasi-particle types in the theory; to identify 
these requires additional physical input (such as their 
self-braiding statistics) . 

It is useful to illustrate how the above counting proce- 
dure works for the case of SU(2) k . Here the coefficients 
are 

i . f{2 3 + iy 



Ci(j) = 



■ sin 



f 



sm 



fc + 2 



simt/(k + 2) \ fc + 2 



(B4) 
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Thus we see that there are indeed two linearly indepen- 
dent combinations of coefficients which will produce the 
desired result. The particle types are given by taking the 
sum and difference of these, to obtain either all even in- 
teger or all odd integer spins on the edges. (For example, 
in SU(2)2, this gives the quasi-particle operators (|, \)\ 
and (|, |)o)- As we have seen by direct computation, 
these are indeed the two expected particle types. 
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